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We show that optically encoded two-qubit Bell states can be unambiguously discriminated with 
a success probability of more than 50% in both single-rail and dual-rail encodings by using active 
linear-optical resources that include Gaussian squeezing operations. These results are in contrast 
to the well-known upper bound of 50% for unambiguous discrimination of dual-rail Bell states 
using passive, static linear optics and arbitrarily many vacuum modes. We present experimentally 
feasible schemes that improve the success probability to 64.3% in dual-rail and to 62.5% in single- 
rail for a uniform random distribution of Bell states. Conceptually, this demonstrates that neither 
interactions that induce nonlinear mode transformations (such as Kerr interactions) nor auxiliary 
entangled photons are required to go beyond the one-half limit. We discuss the optimality of our 
single-rail scheme, and talk about an application of our dual-rail scheme in quantum communication. 

PACS numbers: 03.67.Hk, 42.50.Ex 



Introduction — Bell measurements, jointly projecting 
two qubits onto the so-called Bell basis, constitute a cru- 
cial step in many quantum computation and communi- 
cation protocols, including dense coding [1], quantum 
repeaters [2], and teleportation-based quantum compu- 
tation [3, 4]. The two most common encodings for op- 
tical Bell states are the single-rail (SR) and the dual- 
rail (DR) encodings [5]. An important result with re- 
gards to Bell measurements (BMs) is the impossibility 
of deterministic unambiguous DR BMs using passive lin- 
ear optics [6] , even when arbitrarily many auxiliary pho- 
tons, photon-number-resolving detectors (PNRDs), and 
dynamical (conditionally changing) networks are avail- 
able [7]. Quantitatively, the success probability for un- 
ambiguous DR BMs using a static linear network aided 
by vacuum modes and PNRDs was shown to be tightly 
bounded from above by 50% [8]. In the supplementary 
section of this work, we give a simple proof showing that 
the upper bound of 50% also holds for SR BMs, even 
with the inclusion of dynamical networks [21]. A fully 
deterministic BM requires at least a cubic Hamiltonian 
[22] (which is highly inefficient in practice [9]), unless 
one relies on embedded Bell-state analysis [10]. Current 
proposals for going past the 50% upper bound without 
using experimentally challenging nonlinearities rely on 
using entangled photon ancilla states (which are gener- 
ally expensive and probabilistic to create) and a suffi- 
ciently large interferometer to combine the signal and 
ancilla modes [11], [23]. Similar to Ref. [4], BMs in 
these proposals are made near-deterministic in the limit 
of asymptotically large ancilla states, but without the 
need for conditional dynamics. 

In light of the above facts, the motivation for this 
work is threefold. First, from a theoretical point of view, 
while an upper bound on the success probability of un- 
ambiguous Bell discrimination using passive linear optics 



and auxiliary vacuum modes has been established and 
well celebrated, no such upper bound has been shown 
to exist for active linear-optical circuits that include ad- 
ditional single-mode squeezers and correspond to arbi- 
trary quadratic Hamiltonians [12] (and hence arbitrary 
linear transformations of the mode operators [24]). Even 
though there is a renewed interest in the area of BMs 
[11, 13], there exists a gap in our understanding of what 
is possible between passive linear optics and cubic non- 
linearities [25]. Second, historically, squeezing has typi- 
cally been used in continuous- variable quantum comput- 
ing and information processing [12, 14]. We, however, 
want to explore the use of squeezing in discrete- variable 
BMs, opening up the possibility of combining continuous- 
variable and discrete- variable toolkits for enhanced quan- 
tum information processing [5]. Finally, from an ex- 
perimental point of view, squeezing low-photon-number 
states has become an achievable feat in the last few years 
[15], which warrants an exploration into its potential for 
quantum information processing. Note that in Ref. [15], 
optical squeezing has been promoted from an offline ex- 
perimental resource to a controllable online operation, 
just as is needed for our purposes. 

We show in this paper that beating the 50% bound on 
the success probability of unambiguous BM is possible 
by using single-mode squeezers and beam splitters. Our 
schemes are free of ancilla photons and vacuum modes. 
Hence, the static interferometric networks needed are no 
larger than the signal- mode space (i.e., four modes in DR 
and two modes in SR). We first present a scheme with 
a success probability of 64.3% for DR unambiguous BM 
(while we do not know if this is optimal). We then con- 
sider SR BMs using an arbitrary two-mode network of 
beam splitters and squeezers characterized by real pa- 
rameters. We show that for equal squeezing in the two 
modes, there exists an entire class of interferometers that 
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yield 62.5% success probability. We then look at exper- 
imental considerations and possible applications of our 
schemes. In the supplementary section, we discuss nu- 
merical results that show that it is not possible to beat 
62.5% with a real two-mode active linear network (i.e., 
without phase shifts) for a squeezing of up to 8.686 dB. 

Dual-Rail Bell Discrimination — Let us start with DR 
Bell states in polarization basis, given by 



\iP±) = -L(\HV)±\VH}), 
\4> ± ) = ±={\HH)±\VV)), 



(la) 
(lb) 



where H and V stand for horizontal and vertical polar- 
izations, respectively. Consider the setup shown in Fig. 
[1] with Bell states as the input. We choose a balanced 




FIG. 1: Scheme for DR Bell discrimination that yields a suc- 
cess probability of 64.2%. The Bell states pass through a 
balanced beam splitter (23), two polarizing beam splitters (Vi 
and 7*2)5 an d four squeezers («Si to S4) with a squeezing of 
5.7195 dB each. The output is sent to PNRDs. 

beam splitter whose action on the mode operator vector 



Af- 



4) T is represented by the matrix ^= ^ * 

ter passing through the beam splitter and the polarizing 
beam splitters (where vertically polarized photons are re- 
flected), the Bell states transform into 

|^ + )^^ (11100) + 10011)), (2a) 
|^-)^^ (11010) -10101)), (2b) 

|^) % - (|2000) + |0002) ± |0200) ± |0020)) . (2c) 

We will disregard global phases from now on. Notice 
that at this point the states \^) would be perfectly 
distinguishable if we were to detect four-mode photon- 
number patterns, while the states |0 ± ) would be indistin- 
guishable. Single- mode squeezers add photons in pairs, 
thereby preserving the even and odd parity of the number 
of photons in each output mode. Hence, the states \^) 
remain perfectly distinguishable after applying squeez- 
ing. For this reason, we only concern ourselves with the 
states |0 ± ) in the following analysis. 



Now consider the squeezing operator S given by 
exp [r(a 2 — a) )/2]. After normal ordering, we can write 
its effect on the relevant Fock states |0), |1), and |2) as 



S|0) =Vsech r exp (-tanh ro) /2) |0) , 



<S|1) =(sech r) 6/2 exp (- tanh ra T /2) |1) , 
S\2) ^sech r /2 tanh r exp ( - tanh ra 1 " 2 /2) |0) 



(3a) 
(3b) 



+ (sechr) 5/2 exp (- tanh ra f2 /2) |2) . (3c) 

Each Bell state is transformed into an infinite sum of 
photonic four-mode states when it is passed through the 
squeezers. After some algebra, we can identify the coeffi- 
cients for the different photonic states (the coefficients for 
the states |2000), |4246), etc.). One possible strategy for 
unambiguous discrimination, then, relies on choosing a 
squeezing value where some of these coefficients are zero 
for |</> + ), but non-zero for |</>~), so that some outputs oc- 
cur unambiguously for \4>~). Eq. (4) below shows the 
output from the squeezers up to two-photon terms: 

|^) ->^= |0000) - i (tanh ro^ - sech 4 r)(|2000) + |0002)) 
v2 2 

- -(tanh ra^ T sech 4 r)(| 0200) + |0020)), (4) 

where a ± = tanhr sech 2 r(l ± 1). The complete output 
is presented in the supplementary section, but we can 
already see that the vacuum output is unique to the in- 
put |</> + ), leading to an above-50% success probability 
of BM for any non-zero value of squeezing. By pick- 
ing r = 0.6585 (a squeezing of 5.7195 dB [26]), so that 
tanhra + — sech 4 r = 0, we can make all the two-photon 
output terms unique to \4>~). At this value, the vacuum 
coefficient in |(/> + ) is also at its maximum. Hence, the 
most significant unambiguous output terms for our in- 
put Bell states |^) ± ) (denoted by |^) unique ) are given 
by: 



unique 
unique 



0.5443 |0000) , 

0.2222 |2000) - 0.2222 |0200) 
0.2222 |0020) + 0.2222 |0002) . (5) 



As a by-product of the above condition, a number of 
higher-photon-number terms also become unambiguous 
(i.e., only show up for |</> + )). Adding up the success 
probabilities from these higher-order terms, we obtain a 
total probability of 37.48% and 19.75% for the inputs 
\4> + ) and respectively. Combining this with de- 

terministic discrimination of l^), and assuming equal 
input probability for each of the Bell states, the network 
shown in Fig. [1] yields an overall success probability of 
64.3%. 

Single-Rail Bell Discrimination — In SR, similar to the 
DR case, a scheme relying on a balanced beam split- 
ter followed by two single-mode squeezers, each with a 
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squeezing of 6.2696 dB (r = 0.7218) produces a suc- 
cess probability of 62.5%. However, here we shall con- 
sider a more general setup: a two- mode linear circuit 
parametrized by real reflectivity and squeezing parame- 
ters. The equal-squeezing scheme is later discussed as a 
special case of this general setup. 
Bell states in SR are given by: 



|^> = -J=(|10>±|01», 
|0 ± ) = i=(|OO)±|ll)). 



(6a) 
(6b) 



Using the Bloch-Messiah reduction [16], an arbitrary two- 
mode squeezing network can be decomposed into a com- 
bination of beam splitters and single-mode squeezers, as 
shown in Fig. [2]. The effect of the first beam splitter B\ 




FIG. 2: In SR, we consider the Bloch-Messiah reduction of 
a general two-mode squeezing network. The Bell states are 
passed through a beam splitter Bi(9i), followed by squeezers 
<Si(ri) and £2(72), and finally through a second beam splitter 
82(02); 61 and 62 are the reflectivity parameters, and n and r2 
are the squeezing parameters. The output is sent to PNRDs. 
We only discuss networks with real parameters. 

on the operator vector (a[ a^) 1 is given by the matrix 

cos 0i sin 61 \ ^ ^ e seconc [ beam split- 

- sin &i cos &i J 
ter is given by the substitution 0\ — » 02. Applying the 
beam splitters and squeezers sequentially on the input 
Bell states, the states \^) are transformed into: 



1 . N / sech n sech r 2 , +2 +2 t t 
\W ) ->y 2 exp(xa| + ya\ + za\a' 2 ) 

x (0 ± |10) +a; ± |01>) , 



(7) 



where x = — tanh 7*1 cos 2 02/2 — tanh 7*2 sin 2 02/2, 
y = — tanh ri sin 2 6^2/2 — tanh 7*2 cos 2 #2/2, z = 
cos 02 sin 02 ( — tanh r 1 + tanh r2 ) , and 

= sech 7*1 cos 02 (cos 0i =p sin 0i ) 

— sech 7*2 sin 02 (sin X ± cos 0i), (8a) 
cj ± = sech ri sin 02 (cos #1 =p sin 0\ ) 

+ sech r2 cos 02 (sin #1 ± cos 0i). (8b) 

The case of passive linear optics is obtained by setting 
r = 0. Specifically, |^ + ) -» |01) and \^~) -> |10) can 
be recovered by suitably choosing #1 and #2, e.g., #1 = 
and 02 = 7r/4. By writing out the first few terms in Eq. 



(7), it can be seen that for an apparatus that measures 
photon numbers in each mode, at least one of the four 
conditions (Vt + = 0, Q~ = 0, uj + = or uj~ = 0) must 
be satisfied to unambiguously discriminate between \^) 
with a non-zero probability. If we choose to ignore the 
above constraints, we may distinguish the states \^) 
at most deterministically, but that will yield an overall 
success probability of 50%, which is no better than that 
achieved with a passive linear network. Hence, satisfying 
one of these conditions is necessary to go beyond the 
passive linear upper bound of 50%. 

We now look at the output for the states |^> ± ): 



sech n sech T2 



exp (xa\ + ya\ + za\a\) 



( 7 ±|00)+p ± |02)-p ± |20)+C ± |H)), (9) 



where 



=1 =F sin 2#i (tanh 7*1 — tanhr2)/2, 



(10a) 



± 1/ v 2 (cos 20i sin 202 sech r\ sech T2 
— sin 20i (sin 2 02 sech 2 r\ — cos 2 02 sech 2 7*2)) , 



(10b) 



C = ± cos 20i cos 202 sech n sech T2 



=F sin 20i sin 20 2 (sech 2 7*1 + sech 2 r 2 )/2. (10c) 

Eqs. (7) and (9) show that \^) always result in an odd 
number of photons in the output, while |(/) ± ) result in an 
even number. Hence, we only have to distinguish |^ + ) 
from and |</> + ) from \<j)~). Eqs. (7) and (9) allow 

us to numerically analyze the case of a general two-mode 
network in the supplementary section. In the following, 
we analytically consider the case of equal squeezing (n = 
7*2), for which we achieve a success probability of 62.5%. 

Single-Rail, Equal Squeezing — For 7*1 = 7*2 = r, 
Vt + = => uj~ — 0, and H = => uo + = 0, as can be 
seen from Eq. (8). Hence, imposing any of the four condi- 
tions renders the states \^) completely distinguishable, 
as evident from Eq. (7). Further, we only need to look at 
the constraints Or 1 = to infer the success probabilities 
for uj± = 0. From Eqs. (8) and (10), we have 



ft* = < 



COS 2 2 



(1 ±sin20i)/2 , 
p+ = ±sech 2 r/\/2 , C + = , 
p- = =Fsech 2 r/\/2 , C" = . 



(11) 



Combined with x = y = — tanhr/2, z = 0, and 7* = 1, 
we see that after imposing ft* = the coefficients for 
10*) are independent of 0i and 02- This implies that 
each value of r gives an entire class of networks with 
the same unambiguous-discrimination success probabil- 
ity. For example, a network with a balanced beam split- 
ter after equal-strength squeezers gives the same success 
probability as the one with the order reversed. Further, 
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setting r = gives us an entire class of passive networks 
for which \^) are distinguishable and 10^) are indis- 
tinguishable. As a final observation, since £7 + = and 
Q~ = give only an immaterial sign difference in p ± , 
we are free to choose either condition for determining 
the success probability. In the following, we have chosen 
tt+ = 0. 

We can now evaluate Eq. (9) with the above simplifica- 
tions to write the output |^) ± ) analytically in Fock space. 
The complete output is presented in the supplementary 
section, but the following truncated form highlights the 
main aspects of the calculation (similar to the DR case): 

i , i , sech r . , sech r , , , 9 x . . 

± ) — -=- |00) + (-tanhr Tsech 2 r) |20) 

v2 2 

sech r o 
+ — ^— (-tanhr ± sech 2 r) |02) + (12) 

where we have only kept terms up to two photons. While 
the vacuum output is always ambiguous, the states 1 20) 
and 1 02) can be made unambiguous at r = 0.7218 (a 
squeezing of 6.2696 dB). Adding up the success prob- 
abilities from higher-order unambiguous terms yields a 
success probability of 25.03% for each of the states |0 ± ), 
resulting in an overall probability of 62.5%. We find nu- 
merically that for a phase-free and at most 8.686 dB 
squeezed, but otherwise arbitrary, two-mode squeezing 
network, 62.5% is indeed the best efficiency (discussion 
given in the supplementary part). 

Experimental Success Probabilities — To take into ac- 
count the experimental upper limit on the photon- 
number resolution of the PNRDs, we calculated effective 
success probabilities for equal-squeezing schemes consid- 
ering PNRD resolutions of 2, 5, and 10 photons (while 
still assuming unit detection efficiency). Table [I] shows 
that almost all the benefits of the presented schemes can 
be harnessed with PNRDs that can resolve up to 10 pho- 
tons [27], which may be possible with current technology 
[IT, 18]. 



TABLE I: Experimental Success Probabilities 





PNRD 


Success 


Encoding 


Upper Limit 


Probability 




2 Photons 


26.5% 


Dual-Rail 


5 Photons 


56.1% 




10 Photons 


63.2% 




2 Photons 


41.8% 


Single-Rail 


5 Photons 


56.3% 




10 Photons 


62.0% 



A further important experimental consideration is 
whether, and to what extent, our schemes are robust 
against photon losses (for example, through an inefficient 
detector). The SR scheme, as in the passive case, is not 



robust in this regard. But, note that in DR a passive lin- 
ear setup could detect single-photon losses. Remarkably, 
our squeezing-enhanced DR scheme achieves the same. 
To see this, recall Eqs. (2) and the fact that squeezing 
adds photons in pairs. Hence, after the four single- mode 
squeezers, the Bell states have the form 

|^ + ) — \odd, odd, even, even) + \even, even, odd, odd) , 
\ip~) — ^ \odd, even, odd, even) — \even, odd, even, odd) , 
|0 ± ) — \ev en, even, even, even) , 

where we have only shown the even/odd number of the 
photons in each output mode. The above equation shows 
that measuring an even number of photons in an odd 
number of output modes signals a single-photon loss. 
Such a measurement can be discarded, leading to a re- 
duced success probability depending on the single-photon 
loss rate in the experimental setup. In contrast to the 
passive linear-optical case, two-photon losses cannot be 
reliably identified. 

Applications — Efficient BMs have direct applications 
in quantum computation and communication. Whereas 
measurement-based quantum computation requires at 
least near-deterministic teleportation steps [3, 4], proba- 
bilistic BMs may be used in conjunction with local quan- 
tum memories for connecting the elementary segments 
of a quantum repeater [19]. Such a probabilistic setup 
is sufficient to suppress the exponential decay of entan- 
glement due to channel transmission losses. The pair- 
creation rates in a general repeater chain of total length 
L are proportional to (L/L ) log2 ( 2Pswap / 3 \ with repeater 
stations separated by distances Lo , and the success prob- 
ability for entanglement swapping given by P SW ap- Pswap 
effectively corresponds to the BM efficiency. This rate 
formula is a good approximation for small initial distri- 
bution efficiencies over L and perfect quantum memories 
[19]. For a typical repeater example of L = 5120 km with 
Lo = 20 km, we have L/Lo = 256. This leads to rates of 
~ 0.000152 using optimal, passive BMs, while our local- 
squeezing-enhanced DR scheme (with four squeezers at 
each repeater station) will ideally give ~ 0.001140. Thus, 
in this example, we obtain an improvement of almost 
one order of magnitude, reducing the number of repeater 
chains operating in parallel from 10, 000 to 1, 000 for the 
distribution of at least one long-distance entangled pair 
per unit time. 

Summary and Outlook — We have presented experi- 
mentally feasible schemes relying on static, active lin- 
ear optics that increase the success probability of optical 
qubit BMs beyond 50% without introducing any form of 
ancilla states. Our results open up interesting possibili- 
ties for the use of squeezing in discrete- variable quantum 
information processing. It remains to be seen whether 
the addition of vacuum modes to our schemes increases 
the success probability further, possibly to near 100%. 
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SUPPLEMENTARY MATERIAL 



Throughout the supplement, we use commas to separate the photon numbers in different modes, e.g., what in the 
main text was written as |0002) is written as |0,0,0,2). 



Passive Linear Optics Upper Bound for Unambiguous Bell Measurement in Single-Rail 

The mathematical formulation presented in Ref. [8] (which dealt with dual-rail Bell measurements with the help 
of passive linear optics and vacuum ancillae) could also be applied to the case of single-rail Bell measurements. The 
following argument, however, is not only simpler, but also includes the possibility of conditional dynamics, which was 
excluded in Ref [8]. We show that in single-rail, passive linear networks aided by an arbitrary number of vacuum modes 
and conditional dynamics give a tight upper bound of 50% on unambiguous Bell measurement success probability. 

Consider the input |^) ± ) (g) |£) to a passive linear circuit, where |£) is the (n — 2)-mode vacuum ancilla, |O) 0(n - 2) . For 
concreteness, the states \^) are input in modes one and two (this numbering is arbitrary). The effect of an n-mode 
passive linear circuit is to mix the input creation operators (aj), such that the output operators (b\) are given by 

n 

where Uij is the (i, j)-element of an n x n unitary matrix. This implies that the inputs <g> |£) are mapped to 

n 

|^> ® |0 — (1 ± £ U u a\U 2j a}) \0f n . (14) 

In the above form, it is readily seen that for a static network, the photon-number patterns in the output are identical 
for \4> + ) and \4>~). In fact, the patterns are identical even after allowing for conditional dynamics. To see this, consider 
the cases of detecting zero, one or two photons in some mode s. If we detect two photons in s, then the rest of the 
modes contain the vacuum, which is indistinguishable for the two inputs. Now, say we detect one photon in s. Then 
the conditional state is 

i^ ± )®io-^~±x; a *4io>® (n " i) , as) 

k 

where the summation over k runs from mode one to (s — 1) and from (s + 1) to n, and are some undetermined 
coefficients. Again, the output patterns are identical for the inputs even after further linear processing. To see that 
the case of zero photons in s is also a failure event, notice that after detecting zero photons, the conditional state 
looks similar to that in Eq. (14), except that the conditional state is over n— 1 modes instead of n. This implies, once 
again, that any further linear processing would be futile. Hence, the probability of discriminating between the states 
|0 ± ) is zero. Combined with the fact that a balanced beam splitter deterministically discriminates between the two 
states |V ;± ) 7 we arrive at the tight upper bound of 50% success probability for SR BM using a passive linear-optical 
circuit, an arbitrary number of vacuum modes and conditional dynamics. 



Complete Output in the Dual-Rail Case 



Using Eqs. (3), we can write the effect of the squeezing operator on the states |0) and |2) as 



5|0> =v / sechr ^ 



> y/2n\ ( — tanhr\ n 



n=0 



5|2) 



Vsechr tanh? 

: 71 



n=0 



|2n), 



2n! / — tanhr\ n 



\2n) +sech 5/2 r^ 



tanh r 



n-l 



f2n\ 



V2(n-l)\ 



|2n). 



The next step is to apply the squeezing operator on the states \<j> ) as given in Eqs. (2): 
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, tanhr l2 ^ W '(2m)!(2n)!(2p)!(2g)! /-tanhrY"^^ , 

eech'r J; Vl ^ILi HH (l±l)|2m,2n,2p,2 g ) 



v/2 ^ n m!n!p!<?! 

v m =0,n=0 
p=0,gr=0 

sech 4 r y> y/(2p)l(2q)\ /-tanhr 
2 ^ p\q\ \ 2 J 



, x , , tanhr\ m+n_i v/(2m)!(2n)! , v ^ J(2m)\(2n)\ /-tanhr\ 

x Uul— ) ^^'^^X^UW^i |2 " 2n * 25) 

sech 4 r v /(2m)!(2n)! / -tanhr x m+ ™ 

+ 2 ^ m!n! I 2 

m=0,n=0 v 

/ ~ (-t & xi\ir\ p+q - 1 J(2p)\(2q)\ , 4 ^ /- tanhr \ p+,_1 v / (2p)!(2 ( /)! , ' 
x V V 2m,2n,2p,2g) ± V VI h>) \ hj 2 m,2n,2p,2q) 



Next, we want to identify the output based on the number of photons in each mode. Since this requires forming 
linear combinations of the above summations, the coefficients of the terms turn out to be similar with minor differences 
in relative signs. Hence, in order to write the complete output compactly, we define the functions g ± , h ± , j ± and ^ 
with generic mathematical operations (&i, &2, &3 and h±) as follows: 



sech r y (2m)! / — tanh 



g ± (m J k 1 ) = -^^r 1 { 2 J ( Q &1 m ^ ' (16a) 

^(m,n, &1 , &2 ) = ^V^SS ( ztanhrA^ (q± ^ fe } 
2 m!n! \ 2 / 

±, o o „ x sech 2 r J{2m)\(2n)\{2p)\ /- tanhr \ m+n+p , . c „ „ . , 

j ± (m,n,p ) &i,& 2 ,& 3 ) = — r— — , , o &i m/3 & 2 n/3 & 3 p/3) , (16c) 

fc ±(m,n,p, g , &1 , &2 , &3 ,M = ^ a ^(2m)!(2n)l(2p) !(2g)i / -tanhr^ + ^ ± ^ ^ 

2 m\n\p\q\ \ 2 J v y 

(16d) 

where a ± = \/2 tanhr (1 ± 1) and /3 = —a/2 sech 2 r / tanhr (note that we have defined a ± slightly differently here 
compared to Eq. (4) in order to highlight the exponentially decaying factor of sech 2 r with each coefficient). Using 
these functions we can write the output from the linear network as: 

sech 2 r > 

1^) -> — 2~ « ± |0, 0, 0, 0) + +)(|2m, 0, 0, 0) + |0, 0, 0, 2m)) + 5 ±(m, ±)(|0, 2m, 0, 0) + |0, 0, 2m, 0))] 

m 

oo 

+ [^(m, n, +, ±)(|2m, 2n, 0, 0) + |2m, 0, 2n, 0)) + /^(m, n, +, +) |2m, 0, 0, 2n) + /^(m, n, ±, ±) |0, 2m, 2n, 0) 

m,n 

oo 

+ /i ± (m,n,±,+)(|0,2m,0,2n) + |0,0,2m,2n»] + ^ [^(m, n,p, ±, ±, +) |0, 2m, 2n, 2p) + ^(m, n,p, +, ±, ±) |2m,2n,2p,0) 

m,n,p 

CO 

+ ^(m, n,p, +, ±, +)(|2m, 0, 2n, 2p) + |2m, 2n, 0, 2p))] + ^ A: ± (m, n,p, q, +, ±, ±, +) |2m, 2n, 2p, 2g) , (17) 

m,n,p,q 

where all the lower limits on the summations start from one. Two-photon terms are made unambiguous by set- 
ting their coefficients equal to zero for \<p + ) input, i.e., by setting g + (l,+) = (or equivalently a + + j3 = 0). 
Higher order distinguishable terms can be enumerated by noticing that some of the coefficients h~(m : n, &i, &2), 
j~(m, n : p, &i, &2, &3) and k~(m, n : p, q : &i, &2, &3, &4) in Eqs. (16) and (17) can be made zero for suitable values 
of m, n, p and q. This is succinctly presented in Table [II]: 

Evaluating all the unambiguous terms that give up to 0.01% contribution to the success probability, we arrive at a 
probability of 37.49% for |</> + ) and 19.75% for \4>~) 1 giving us an overall success probability of 64.3%. 
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TABLE II: Higher Order Unambiguous Terms in 



Condition Unambiguous Terms in |</> + ) 



m 


— n 


|2ra,2n,0,0) , 


|2ra,0,2n,0) , |0, 2m, 0, 2n) , |0,0,2ra,2n) 


m 


— n + p = 


|2ra,0,2n,2p) , 


|2ra,2n,0,2p) 


m 


— n — p = 


|2ra,2n,2p,0) 




m 


— n — p + q — 


|2ra,2n,2p,2<?) 





Complete Output in the Single-Rail, Equal Squeezing Case 

For the case of equal squeezing, the complete \<tr^) output can be written as: 



tanhr =p nsech 2 r) |2n, 0) + (— tanhr ± nsech 2 r) |0, 2n)) 



tanh r ±m sech 2 r =p n sech 2 r) |2n, 2m) j , (18) 

where the lower limit on the summations again starts from one. To set the coefficient of the 1 02) term equal to zero 
for input, we must satisfy the condition — tanhr + sech 2 r = 0, which implies r = 0.7218 (a squeezing of 6.2696 
dB). This condition is also fulfilled for higher photon-number states of the form |2n, 2m) whenever m — n — ±1. For 
example, the states |2,4) and |4, 6) unambiguously distinguish \<j) + ) with a probability of 1.2% and 0.1%, respectively. 
Similarly, by symmetry of the above equation, the states |4, 2) and |6,4) unambiguously distinguish \4>~) with a 
probability of 1.2% and 0.1%, respectively. Adding up the contributions from all the terms that give at least 0.01% 
success probability gives us an overall probability of 62.5%. 



i«^(io,o> + E(^rr^((- 

sechr /y, / -tanhr \ m+ "" 1 V(2m)!(2n)! ,_ 
+ \/2 \^\ 2 ) 2m\n\ ^ 

v \m,n x 7 



Numerical Calculations for the Single-Rail Case 

Numerical calculations on a two-mode network with real parameters were done in Mathematica (the code is available 
upon request). Quantum Mathematica add-on by Jose Luis Gomez-Munoz and Francisco Delgado was used for part 
of the simulation, which is available at http://homepage.cem.itesm.mx/lgomez/quantum/. The objective was to 
see if the success probability of 62.5% in SR is optimal within the context of a two-mode network characterized by 
real parameters in the B loch- Messiah Reduction. 

Eqs. (7) and (9) were used to numerically evaluate the success probability. We swept 6i in steps of 0.017T radians, 
and the unit-less parameters r\ and r<2 in steps of 0.05 from to 1 (for a maximum squeezing of 8.686 dB). $2 was 
constrained by one of the four conditions £7 + = 0, Q~ = 0, uj + = or uo~ = 0. Bell states were simulated to have up 
to 26 photons in each mode (implying an error of less than 1% from neglected terms in the success probability at a 
squeezing of 8.6859 dB). An output that was 400 times more likely to occur for an input state \A) than for any other 
input was considered to unambiguously identify \A). 

Calculating the success probability for a total of 352, 800 numerical data points, it was found that, within a numerical 
accuracy of 1%, the optimal success probability was 62.5%, which happens for equal squeezing in the two arms of the 
network between 6.08 dB and 6.51 dB. This is shown in Fig. [3] where we plot the success probability as a function 
of squeezing and 6\ for equal squeezing in the two arms (#2 stays constrained by the condition Vt + = 0). 

In Fig. [4] , we present success probability plots as a function of squeezing in the two arms of the linear network for 
randomly chosen values of 6\ = 0, 0.677T, 1.357T and 1.87T radians. These plots are all for the data set constrained by 
the condition ^ + = (the other three conditions produce similar plots). We would like to point out that the plots 
are not symmetric in the squeezing of the modes, since the coefficients given in Eq. (10) are not symmetric in r\ 
and T2- As can be seen from the plots, qualitatively, the peak value of the success probability is about 60%. This 
can be quantitatively corroborated by looking at the success probability of individual data points. Since the plots 
do not add much to an intuitive understanding of the interplay between active and passive linear optics, we did not 
consider it worthwhile adding plots for other values of 9\. It is unlikely that generalizing the numerical calculations 
to include phase shifts in an arbitrary two-mode linear network will change the upper bound, though this will have 
to be addressed in future work. 
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FIG. 3: Success probability (P) for n = T2- #2 is constrained by Q + =0. P is bounded from below by 50% since the states 
|0 ± ) are deterministically distinguishable for n = r^. Since #2 is constrained, every value of 61 gives a success probability of 
62.5% at a squeezing of 6.2696 dB. 




(a) 0i = radians (b) 6\ = 0.677T radians 




(c) 61 = 1.357T radians (d) 6\ = l.Sir radians 



FIG. 4: Success probability as a function of n and T2 for various values of Q\ after imposing the condition Q + = 0. Si and 
S2 stand for squeezing in modes one and two of the two-mode network, respectively. For every value of 61 we see a peak in 
the discrimination success probability for equal squeezing in the two modes of around 6 dB, as expected from the analytical 
calculation of the equal squeezing case presented in the main text. 



